Ultraslow Propagation of Matched Pulses by Four- Wave Mixing in an Atomic Vapor 
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We have observed the ultraslow propagation of matched pulses in nondegenerate four- wave mixing 
in a hot atomic vapor. Probe pulses as short as 70 ns can be delayed by a tunable time of up to 
40 ns with little broadening or distortion. During the propagation, a probe pulse is amplified and 
generates a conjugate pulse which is faster and separates from the probe pulse before getting locked 
to it at a fixed delay. The precise timing of this process allows us to determine the key coefficients 
of the susceptibility tensor. The presence of gain in this system makes this system very interesting 
in the context of all-optical information processing. 

PACS numbers: 42.50.Gy, 42.65.Yj 



Slow group velocities, valuable for all-optical signal 
processing, are obtained at a resonance peak of the trans- 
mission spectrum of a medium, and a number of differ- 
ent implementations of this principle have been demon- 
strated. They rely either on a reduction of the absorp- 
tion, such as electromagnetically induced transparency 
(EIT) coherent population oscillations Q, and dual 
absorption lines or on a gain resonance, like stimu- 
lated Brillouin scattering [3] and stimulated Raman scat- 
tering [5] . To be useful in the context of all-optical signal 
processing, an optical delay line should be able to pro- 
duce a fractional delay (defined as the ratio of the delay 
to the duration of the pulse) larger than unity with only 
modest absorption and pulse broadening. Recent devel- 
opments 0, H have shown the benefits of using an 
amplifying medium to alleviate the absorption and dis- 
tortion issues usually associated with slow light Q. 

We have examined the group velocity reduction effects 
due to nondegenerate four-wave mixing (4WM) in hot ru- 
bidium vapor, and have obtained large fractional delays 
with almost no distortion. The presence of gain in this 
system makes it in principle possible to stack such delay 
lines and achieve fractional delays only limited by pulse 
broadening. Another notable feature of the amplifica- 
tion in the 4WM process is the generation of a conjugate 
pulse which is coupled to the probe and which propagates 
alongside it, similar to matched pulses in EIT systems Q. 
We have studied the interplay between the 4WM coupling 
of the probe and conjugate, and the Raman coupling of 
the probe and pump. This interplay leads to the ultra- 
slow propagation of matched probe and conjugate pulses 
and to the enhancement of the 4WM gain. Such an en- 
hancement, when applied to cross-phase modulation, is 
the key element of recent optical quantum information 
processing proposals 1(J 11 , 
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Our apparatus, which is essentially the same as the 
one described in Ref. consists of a linearly polar- 
ized, continuous, strong (up to 280 mW) pump and a 
cross-polarized, pulsed, weak (0.5 mW) probe propagat- 
ing at a small angle (0.5°) through a 2.5 cm-long 85 Rb 
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FIG. 1: (color online), (a) Experimental setup. PBS: polar- 
izing beamsplitter, (b) Energy-level diagram of the DI line 
of 85 Rb, showing the double-lambda scheme. Note that the 
pumps fii and 0,2 are in fact the same laser beam, (c) Probe 
transmission profile versus two-photon detuning S. 



cell heated to 90°C-140°C (Fig.QJi). The pump and the 
probe are near resonant with a Raman transition between 
the two hyperfine electronic ground states of 85 Rb, with 
a controllable detuning S (Fig. [TJd), and have 1/e 2 radii 
of 600 /im and 350 fim, respectively. The residual 2- 
photon Doppler broadening due to the small angle is a 
few MHz. The detuning from the 5Pi/2 excited state is 
Ai/27r w 850 MHz, and the peak pump intensity (up 
to 45 W/cm 2 ) is high enough to excite off-resonant Ra- 
man transitions with a detuning A/2-7T rs 4 GHz from 
the excited state. The double-lambda [3| is closed by 
the conjugate beam which emerges on the other side of 
the pump from the probe, with the same polarization 
as the probe. The combination of the beam polariza- 
tions and the Zeeman substructure makes the system a 
four-level system (the two virtual excited states are or- 
thogonal). The probe amplification is sharply resonant 
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in S, as shown in Fig. QJ:, with a gain that can reach 30. 
The gain feature leads to a strong dispersion of the index 
of refraction and thus a low group velocity for the probe. 

We measure the group velocity delay by recording the 
arrival time of a 70 ns-long (full width half maximum 
[FWHM]) gaussian probe pulse with and without the 
atomic medium (reference pulse). Figured shows an ex- 
ample in which the parameters are set to provide a large 
probe gain (G = 13). By varying the two-photon detun- 
ing 8 and the pump intensity, one can tune the probe 
delay and achieve a fractional delay larger than 0.5 such 
that the pulse remains gaussian and is broadened by less 
than 10% of its original width. This low level of distor- 
tion is remarkable in comparison with that seen in some 
EIT experiments 0] . The maximum delay corresponds to 
a group velocity of c/500, where c is the speed of light in 
vacuum. It is in general possible to tune the pump inten- 
sity, the pump and probe detunings, and the temperature 
to achieve an overall gain of unity. Figure [2k also shows 
the record of the conjugate intensity. A striking feature 
is the emergence of the conjugate pulse before the probe 
pulse. This relative delay is a fundamental property of 
the dynamics of the system, and was predicted and ob- 
served in Refs. [HI, Ha] in the case of resonance on the 
"lower" lambda transition (Ai = 0). 
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FIG. 2: (color online), (a) Slow-light effect near the peak 
of the resonance. The reference pulse is magnified 13 times. 
Thin black solid and dashed lines: probe pulses for the param- 
eters (8/2n, pump power) equal to (10 MHz, 280 mW) and 
(22 MHz, 200 mW) respectively. Color lines: corresponding 
conjugate pulses. The pulses (probe and matching conjugate) 
are broadened by 5% in the less retarded case, and 10% in the 
more retarded case, (b) For a detuning Ai closer to the Ra- 
man absorption dip in Fig.QJ, large delays and pulse breakups 
are observed (here for the conjugate). 

Larger delays, shown in Fig. [^Ja, can be achieved by 
setting 6 between the gain peak and the Raman absorp- 
tion dip present at 6 < (see Fig. [Th). The competition 
between large amplification and large absorption leads 



to some complex dynamics which can result in pulse 
breakup, in a similar fashion to the dual-field solitons 
predicted to exist in three- level systems [17[. In this pa- 
per, we restrict ourselves to the regions of low absorption 
where, in spite of complications associated with the use 
of hot atoms, our system can be consistently described 
over a broad range of parameters by the simple concepts 
developed in the theory of Refs. [3, [lfl . 

We neglect the hyperfine splitting of the excited state. 
Averaged over the Zeeman substructure, the dipole mo- 
ments of all four transitions are equal, which gives fl\ 
and f^2, the peak resonant Rabi frequencies of the pump 
for the "lower" and the "upper" lambda respectively, the 
same value, denoted Q. The double-lambda system in its 
ideal incarnation operates in the limit Ai <C A and has 
the following crucial features. First, a ground state coher- 
ence is established by the "lower" , more resonant lambda. 
The coherence has a lifetime l/j c , limited by magnetic 
fields, collisions and the transit time in the laser beams, 
and corresponds to a dark state in which the absorp- 
tion of the probe is reduced (EIT). Second, the "upper", 
less resonant lambda slightly perturbs this coherence and 
creates a resonant atomic polarization at the probe and 
conjugate frequencies via 4WM, while keeping the pop- 
ulation in the excited state near zero. The dynamics of 
the system can thus be broken down into two intertwined 
processes: EIT and 4WM. 

In the limit of a strong pump and low pump deple- 
tion, most of the atomic population is in the ground state 
F = 3, and the Fourier components £ p (u>) and £ c (— us) 
of the slowly-varying envelopes of the probe and conju- 
gate fields (of wavevectors k p and k c ) obey the equa- 
tions [H, [III : 



(iuj + cd z )E p = ii]A R £* — rj 
(iw + cd z )£* = -ir]A R £ p . 



7c 



(2) 



We assume perfect phase matching. In the limit of 
fi 2 /4Ai > 5, 7c, and A > Ai, 7 (where 7/2-7T = 6 MHz 
is the linewidth of the atomic transition), the coefficients 
in Eqs. © and © are [3 77 = g 2 N/[fl 2 /A + Ai(tf + w + 
he)} ~ 4g 2 N/n 2 = c/v g > 1 and A R = f2 2 /4A. Here, 
g 2 = ckp 2 /(2eoh), k = k p w k c , N is the atomic density, 
and p is the average dipole moment acting on the probe 
and the conjugate. 

The interpretation of Eqs. ([T]) and @ is straightfor- 
ward. The probe field £ p is slowed down by a factor 77 
due to the EIT interaction with the pump [second drive 
term on the right-hand side of Eq. ([TJ]. The EIT reso- 
nance is light-shifted by the pump on the "upper" lambda 
and occurs at S = S — = 0. In addition, £ p and £ c 
are cross-coupled with a coupling constant a — rjAji, re- 
sponsible for the 4WM amplification. The presence of 
X] highlights the role of the longer interaction time due 
to the slow-down effect in obtaining a sizeable nonlinear 
coupling ll|. The other factor in a, the so-called Ra- 
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man bandwidth A#, is the Rabi frequency of a fictitious 
resonant Raman transition driven on both legs by the 
pump field and with an intermediate Raman detuning 
A. As shown by the absence of any dependence on A 1; 
the 4WM dynamics is dominated by the "upper" lambda, 
which acts as a bottleneck. The propagation equations 
are asymmetrical. The imaginary part of the coefficient 
of the direct term for the conjugate [Eq. (j^J], correspond- 
ing to a slow-light effect, is negligible compared to the 
same term for the probe because A ^> Ai 18]. The 
real part, corresponding to a Raman amplification, scales 
as 1/A 2 and is negligible compared to the cross-term a 
which scales as 1/A [l6l. fl8|. As a result, in the absence 
of the 4WM coupling ("bare" fields), the probe and the 
conjugate propagate at velocities v g and c, respectively. 
The finite decoherence j c translates into a small absorp- 
tion of the probe. 

Our system departs from the ideal case described by 
the expressions of rj and a given above in many resp ects. 
Unlike the experiments described in Refs. [Tig, EH, EH, 
which were performed with a resonant probe and a weak 
pump, our probe is tuned to the side of the Doppler pro- 
file, in an already almost transparent region. As a result, 
the position of the gain peak is not tied as closely to 
a narrow EIT window. It depends on the balance be- 
tween the losses, which include the Raman absorption 
dip and the absorption from the Doppler broadened 1- 
photon transition, and the 4WM gain. Factors influenc- 
ing the peak position are the spread of values for Ai due 
to the Doppler broadening, the spread of values for Sli 
and f?2 due to the Zeeman degeneracy, the contribution 
of the usual dispersion of the Doppler broadened vapor 
to the slow-down of the probe, and the only approximate 
phase matching. In practice, this means that the posi- 
tion of the gain peak varies by up to 20 MHz depending 
on parameters like the temperature, Ai and the probe 
intensity. 

In spite of the added complexity and the difficulty of 
directly calculating 77 and Ar, we assume that the prop- 
agation equations (p} and ((2]) are still valid over most of 
the resonance peak, provided that the peak is at S ps 0, 
and that j c <C 2A#. Solving them in the limit 77 ^> 1 
and starting from a probe field £o(&) and no conjugate 
field, one finds: 



£ p (u},z) = £o(uj) exp (i<j(uj)— ^ 
cosh U(u>)- ) 



-M-f^sum £H- 
cj f w V c 



(3) 



= £ Hexp( i aM^)|^sinh(eH^)(4) 

where £ (uj) = ^ a(uo) 2 — a(uj) 2 and a(uj) = , -^-(S + uj + 
ijc). Since j c <C 2Ar, £(uj) is real. As expected, past 
an initial linear growth of the conjugate, both fields grow 
exponentially in distance with (£ — ^7 c )/c as the linear 



gain coefficient. In the limit 5 = 0, one has £ w a. 

The solutions also contain information about the prop- 
agation delay. Equations © and ^ show that the 
fields accumulate a phase across the medium, denoted 
9{u). The main contribution to the group delay -^d{uj) 
for both fields comes from the first exponential and 
gives a common delay r = ^Re[cr(tj)].z/c| ( j = 
rjz/2c. In other words, the probe and the conjugate 
are slowed down by half the "bare" slow-down factor 
7/. The probe experiences an extra delay due to the 
second term in Eq. ([3]). At large gain, the cosh and 
sinh functions are equal and the additional delay is 
At = £Re[a( W )/£( W )]/(l-Im[a(c)/£(c)])| s=0 ^ =0 = 
r\l [2£ — 777J ~ 77 / 2^. At low gain, a first order expansion 
in z gives At = t]z/2c. The picture emerging from this 
analysis is the following: the conjugate pulse is created 
without delay by the probe and travels at a velocity 2v g 
(<C c). The probe pulse travels initially at a velocity v g 
and then locks onto the conjugate by accelerating to 2v g 
when the delay reaches rj/2^ (at a gain close to 2). 
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FIG. 3: (color online). Two-photon detuning scan at a tem- 
perature of 140° C, using 120 ns-long pulses. The bare state 
2-photon resonance corresponds to 8 = 0. The calculated 
gain, inferred from the delays, assumes linear losses equal to 
14% of the peak linear gain. 

We test this interpretation by first scanning the two- 
photon detuning 6, using 120 ns-long (FWHM) probe 
pulses and a pump power of 280 mW, which corre- 
sponds to a spread of ui of 10 MHz around zero and to 
n/2ir = 420 MHz. Figure [3] shows the measured gain, 
the probe delay r + At, which includes the contribution 
of all the retardation effects, and the differential delay 
At. The contribution to the probe delay of the usual 
dispersion effect, which is measured with the pump in- 
tensity strongly reduced, is found to be 8 ± 2 ns. From 
the experimental data and the theoretical expressions of 
r and At, one can deduce a value of 77 and £ for each S, 
in the limit of large gain and small 7 C . Inserting these 
values into Eq. §5§ (with 6 w 0) and adjusting 7 C to O.57 
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(making the linear losses equal to 14% of the peak lin- 
ear gain) 22] , one can reproduce the gain curve with 
reasonable accuracy. For our parameters, the EIT reso- 
nance is light-shifted to 6 — 11 MHzx27r rj 2j, close to 
the observed gain maximum (5 = 6 MHzx27r ~ 7). It 
can be shown that in the above calculations, the approx- 
imation 6 is valid as long as 5 -C 2A# (47 for our 
beam parameters). For \S\ larger than a few 7, that is to 
say in the wings of the gain peak, the approximation is 
expected to break down. 

Next, we directly observe the locking between the 
probe and the conjugate during the propagation. It is 
impractical to continuously vary the distance of propa- 
gation and we instead vary the atomic density N via the 
temperature, which is equivalent. Indeed, a and £ are 
proportional to N through their dependence on r), and 
changing N is like renormalizing z in the solutions ([3]) 
and ((4}. According to ((3|), the renormalized propaga- 
tion length L is related to the probe intensity gain G by 
L = cosh _1 (-\/G). The detuning 6 is set to 15 MHzx2tt, 
near the gain maximum, the pump power is still set to 
280 mW, and the measured delays as a function of the 
renormalized distance for a temperature scan of 50° C 
around 120°C are shown in Fig. [4] The two regimes of 
propagation are very clear. First, the pulses separate in 
time, and second they lock to each other at a fixed de- 
lay. The time offset at the origin is not well understood. 
A direct evaluation of 2 A/ft 2 w At in the ideal case 
using our beam parameters gives a value of 7 ns, compa- 
rable to the one measured near the gain maximum (see 
Fig. [3]). We checked qualitatively that At increases when 
the pump intensity decreases. It is worth noting that the 
detail of the low-gain transient regime depends on the 
initial conditions. For instance, swaping the frequencies 
of the probe and the conjugate would lead to an initial 
propagation in which the probe travels at a velocity c 
while the conjugate travels at a velocity 2v g . 

Finally, an important feature of the model is that the 
gain saturates with the pump intensity. For a gain peak 
location 6 and a decoherence 7 C of the order of 7, £ sat- 
urates when f2 ^> 2^/A-f = 300 MHz x2-7r. In agreement 
with this prediction, we observe that G starts saturating 
at our operating intensity. 

To conclude, we have observed the ultraslow propaga- 
tion of probe and conjugate pulses with matched shapes 
and group velocities through a rubidium vapor. The 
study of the coupled propagation gives access to the 
atomic dynamics through a simple model that reflects a 
few key concepts. Although the hypothesis of the model 
does not match precisely the conditions of the experi- 
ment, our findings on slow propagation and delay lock- 
ing of the probe and conjugate pulses are generic to the 
double-lambda system. The quality of the retardation 
effect in terms of fractional delay and absence of loss and 
distortion suggests the possible existence of a dual-field 
soliton, which would be the result of higher order terms 




FIG. 4: (color online). Probe delay t + At, conjugate delay 
t, and differential delay At as a function of a pseudo propa- 
gation distance. 



in the propagation equations. Such a soliton has been 
predicted in related 4WM schemes [l||. 

It should also be noted that this double-lambda scheme 
is known to generate relative-number squeezed twin 
beams 13] (as well as correlated photons [20J). Extend- 
ing this semi-classical pulse theory to the quantum cor- 
related beams raises two comments. First, as pointed 
out in Ref. [lfij . the time lag between the probe and the 
conjugate should be the limiting factor to the squeez- 
ing bandwidth observed in Ref. [13|. Second, the system 
could be used with gain close to unity to slow light in the 
quantum regime, while preserving nonclassical correla- 
tions, possibly more efficiently than with EIT alone 21 1 . 
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